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Abstract. An important problem in the theory of finite dynamical systems is to link the structure 
of a system with its dynamics. This paper contains such a link for a family of nonlinear systems 
over an arbitrary finite field. For systems that can be described by monomials, one can obtain 
information about the limit cycle structure from the structure of the monomials. In particular, the 
paper contains a sufficient condition for a monomial system to have only fixed points as limit cycles. 
The condition is derived by reducing the problem to the study of a Boolean monomial system and 
a linear system over a finite ring. 



1. Introduction 

Finite dynamical systems are time-discrete dynamical systems on finite state sets. Well-known 
examples include cellular automata and Boolean networks, which have found broad ap plic ations in 
engineering, computer science, and, more recently, computational biology (See, e.g., BMU 
for biol ogic al applications.) More general multi-state systems have been used in control theory 
[ill : Eol : E3 : E3 |. the design and analysis of computer simulations 0; El S 03- One underlying 
mathematical question that is common to many of these applications is how to analyze the dynamics 
of the models without actually enumerating all state transitions, since enumeration has exponential 
complexity in the number of model variables. The present paper is a contribution toward an answer 
to this question. 

For our purposes, a finite dynamical system is a function / : X — * X, where X is a finite set 
[l7) | . The dynamics of / is generated by iteration of / and is encoded in its phase space S(f), 
which is a directed graph defined as follows. The vertices of S(f) are the elements of X. There is 
a directed edge a — > b in S(f) if f(a) = b. In particular, a directed edge from a vertex to itself is 
admissible. That is, S(f) encodes all state transitions of /, and has the property that every vertex 
has out-degree exactly equal to 1. Each connected graph component of S(f) consists of a directed 
cycle, a so-called limit cycle, with a directed tree attached to each vertex in the cycle, consisting 
of the so-called transients. 

Any Boolean network can be viewed as a finite dynamical system / : F2 — > F2, where F2 
is the finite field on two elements and n > 1. In this paper, we study finite dynamical systems 
/ : Fg — ► F™, where F g is any finite field and |F ? | = q. To be precise, we present a family of 
nonlinear finite systems for which the above question can be answered, that is, for which one can 
obtain information about the dynamics from the structure of the function. 

Let / : F™ — > F™, n > 1 be a finite dynamical system. Observe that / can be described in 
terms of its coordinate functions fa : F™ — ► W q , that is, / = . . . , /„). It is well known that any 
set-theoretic function g : W q — > W q can be represented by a polynomial in W q [xi,... ,x n ], see |2l|, 
pp. 369]. This polynomial can be chosen uniquely so that any variable in it appears to a degree 
less than q. That is, for any g there is a unique h € F 9 [xi, . . . , x n ]/(xf — X{\% = 1, . . . , n), such that 
g(a) = h(a) for all a £ F™. Consequently, any finite dynamical system over a finite field can be 
represented as a polynomial system. This is the point of view we take in this paper. 
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In the case that all the fi are linear polynomials without constant term, the dynamics of the 
linear system f can completely be determined from its matrix representation [121 ]. Let A be a 
matrix representation of a linear system / : F™ — ► F™. Then the number of limit cycles and their 
length, as well as the structure of the transients, can be determined from the factorization of the 
characteristic polynomial of the matrix A. The structure of the limit cycles had been determined 
earlier by Elspas [lOj], and for affine systems by Milligan and Wilson |2J]. 

In this paper we focus on the class of nonlinear systems described by special types of polynomials, 
namely monomials. That is, we consider systems / = {fi), so that each fi is a polynomial of the 
form x^x®2 • • • x H r i or a constant. Without loss of generality we can assume that no coordinate 
function is constant, since the general case is easily reduced to this. Some classes of monomial 
systems and their dynamic behavior have been studied before: Monomial cellular automata |5cll4l|. 
Boolean monomial systems IjJ, monomial systems over the p-adic numbers [16c l25| , and monomial 
systems over finite fields (1; 0] . 

In | a special class of Bodean monomial systems was studied, namely those which have only 
fixed points as limit cycles, so-called fixed point systems. The motivation for considering this class is 
the use of polynomial systems as models for biochemical networks. Depending on the experimental 
system considered, such networks often exhibit steady state dynamics. That is, their dynamic 
models have phase spaces whose limit cycles are fixed points. For the purpose of model selection 
it would be useful to have a structural criterion to recognize fixed point systems. The main result 
of the present paper is to reduce this question for monomial systems over a general finite field F g 
to the same question for an associated Boolean monomial system and a linear system over a ring 
of the form 1i/{q — 1). 

2. Reduction of monomial systems 

Let / = (/i, . . . , f n ) : Fg — ► F™ be a polynomial system, with each fi a monomial, that is, 
fi = x^x^' 2 • • • ij", with a%j a nonnegative integer. That is, / can be described by the exponent 
matrix A = (a y -). We first associate with / a Boolean monomial system T{f) and a linear system 
L{f) over the ring R = Z/(g — 1). Recall from [^] that / is called a fixed point system if all limit 
cycles of / consist of a fixed point. We will show that / is a fixed point system if and only if T{f) 
and L{f) are fixed point systems. 

Definition 2.1. For u = {u\ , . . . , ?%) G F™, we define support u, denoted by supp(u), to be 
v = {vi,. . . ,v n ), where 

'1 if Ui ^ 0; 
if Ui = 0. 

The monomial system f = {fx, ■ ■ ■ , fn) induces a Boolean monomial system T{f) = {g\,..., g n ) 
on Fg by setting gi = x\ A ■ ■ ■ x% n , where fi = x" 11 • • • x™ in and v = supp(u). 

Lemma 2.2. There is a commutative diagram 

i i 

Proof. The proof is a straightforward verification, since 

supp(/(u)) = /(supp(u)). 

□ 



Since / = T{f) on the set of all u such that supp(u) = u, the following corollaries are straight- 
forward. 
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Corollary 2.3. The phase space ofT(f) is a subgraph of the phase space of f. 

Corollary 2.4. Suppose that T(f ) is a fixed-point system. If {u, /(u), . . . , /*(u) = u} is a cycle 
in the phase space of f , then supp(u) = supp(/ l (u)) for all 1 < i < t. 

For more results in this direction, see 

Example 2.5. Let 



/ = (xlx 2 ,x 2 xl,xlx 2 X3) : Fg 



Figures n and HI display the phase spaces of the system / and its " Booleanization" T(f), respec- 
tively. 




Figure 1. The phase space of /. 




Figure 2. The phase space of T(f). 

Next we associate to / an n-dimensional linear system over a finite ring. Observe first that 
F* = W q — {0} is isomorphic, as an Abelian group, to Z/(q — 1) via an isomorphism 

log : F* — » Z/(q - 1), 

given by the choice of a generator for the cyclic group F* . Note first that the set of vectors u G F™ 
with all non-zero entries is invariant under /. 

Let a be a generator for the cyclic group F* = Z/(q — 1), and let 

u = (ui,...,u n ) = (a sl ,...,a Sn ) = a s . 

Then /(u) = f{a s ) = a A s . 

Definition 2.6. Define L(f) : (Z/(g - l)) n — ► (Z/(g - l)) n by 

L(f)(a s ) = a A - s . 

As defined, L{f) is a linear transformation of Z-modules. But we can consider it as a linear 
transformation of Z/(q — l)-modules, viewing Z/(q — 1) as a (finite) ring, which we denote by R. 
That is, we have the linear transformation 

L(f) = A: R n — ► R n . 

The proof of the following lemma is a straightforward verification. 



4 O. COL6N-REYES, a. s. jarrah, r. laubenbacher, and b. sturmfels 



Lemma 2.7. There is a commutative diagram 

( F *)« 1+ (F *)n 
I I 

Note that the vertical arrows are isomorphisms. This implies that they preserve the phase space 
structure, including the length of limit cycles JjJ. In particular, we have the following corollary. 

Corollary 2.8. The phase space of L{f) is isomorphic to the subgraph of the phase space of f 
consisting of all states with support vector (1, 1, . . . , 1). 

Example 12.51 (Cont.). For the monomial system / in Example 12. 5| L(f) : (Z/2) 3 — ► (Z/2) 3 is 
defined by L(/)(s) = A • s, where 

10 
.1-1 1 
1 1 



The phase space of L(f) is given in Figure El 




Figure 3. The phase space of L(f). 

We now come to the main result of this section. 

Theorem 2.9. Let f : F™ — > F™ be a monomial dynamical system. Then f is a fixed point system 
if and only ifT(f) and L(f) are fixed point systems. 

Proof. From Corollaries 12.31 and 12.81 if / is a fixed point system, then so are T(f) and L(f). For 
the converse, we assume that L(f) and T{f) are fixed point systems, but / is not. For each limit 
cycle of /, either all states involved have all coordinates non-zero, or all states involved have at 
least one zero coordinate. In the first case it follows that L(f) has a limit cycle of the same length. 
Hence, if / has a limit cycle of length greater than 1, then it can involve only states with at least 
one zero coordinate. 

Let ai, ... , at be the states in the limit cycle. Since this limit cycle has to project to a fixed 
point of T(f) it follows that the a^ have the same support vector, i.e., the same pattern of zero 
entries, and differ only in the non-zero coordinates. Furthermore, the monomials in the nonzero 
coordinates do not involve any variables corresponding to zero coordinates. Thus, if we construct 
new states bj by replacing each in a^ by a 1, the bj will be part of a limit cycle of length at least 
t, which is a contradiction. This completes the proof of the theorem. □ 

3. Linear Systems over Finite Commutative Rings 

The theorem in the previous section shows that to decide whether a given monomial system 
/, over a finite field F ? , is a fixed point system it is sufficient to decide this question for an 
associated Boolean system, for which a criterion has been developed in [gj], and a certain linear 
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system over a finite ring Z/m. It therefore remains to develop a criterion for linear systems over 
finite commutative rings that helps decide whether the system is a fixed point system. Here we 
reduce the case of a general Z/m to that of m being a prime power. 

Let m = st for relatively prime integers s and t, and let / be a linear system over Z/m of 
dimension n. Choosing an isomorphism Z/m = Z/sxZ/(we see that / is isomorphic to a product 

g x h : (Z/s)™ x (Z/t) n — » (Z/s) n x (Z/t) n , 

where 5 and h are linear systems over Z/s and Z/i, respectively. Using the fact that the phase 
space off is then the direct product, as directed graphs, of the phase spaces of g and h (see, e.g., 
[L2| or |l3j), we obtain the following result. 

Proposition 3.1. Let m = st for relatively prime integers s and t, and let f be a linear system 
over Z/m of dimension n. Let g and h be the induced linear transformations over Z/s and Z/t, 
respectively. Then f is a fixed point system if and only if g and h are fixed point systems. 

For the purpose of developing a criterion to recognize fixed point systems it is therefore sufficient 
to study linear systems over rings of the form Z/p r for primes p. The following theorem provides 
a criterion for a further reduction of the problem to a linear system over the prime field Z/p. 

Theorem 3.2. Let f = (fi,...,f n ) '■ (Z/p r ) n — ► (Z/p r ) n be a linear map, and let g be the 
projection map of f on Z/p. That is g = (<?i, . . . , g n ) : (Z/p) n — ► (Z/p) n , where gi = fi mod p. 
Then the phase space of g is isomorphic to a subgraph of the phase space of f . 

-1 



Proof. Let a : {Z/p) 1 



(Z/p r ) n be given by a (a) = a(a±, . 



{CLIP 



,a n p 



- 1 ) 



ap r . Then it is easy to check that a o g = f o a, since the fi are linear maps for all i. Therefore, 



it is straightforward to check that g(a) = b if and only if f(ap r 
space of g is isomorphic to subgraph of the phase space of /. 



hp r j and hence the phase 

□ 



Corollary 3.3. Let f and g be as above. If g is not a fixed point system, then f is not a fixed 
point system. 

The dynamics of projection maps has been studied in jq|. 

Example 3.4. Let / : (Z/4) 2 — ► (Z/4) 2 be given by / = (2a?i + 3x 2 ,x 1 ). Then g = (x 2 ,x 1 ) : 
(Z/2) 2 — > (Z/2) 2 . The phase spaces of / and g are in Figures |U and 03 respectively. 






Figure 4. The phase space of /. 
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Figure 5. The phase space of g. 

It remains to study the dynamics of linear systems over finite rings, in particular to find a 
criterion for a linear system to be a fixed point system. Together with the results in this paper, 
such a criterion would provide an algorithm for deciding whether a monomial system over an 
arbitrary finite field is a fixed point system. However, it appears to be a difficult problem to 
understand the dynamics of linear systems even over rings of the form Z/p r , due to the lack of 
unique factorization in the polynomial ring Ti/p r [x}. See, e.g., 
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